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We reanalyze transfer matrix and Monte Carlo results for the critical Binder cumulant U* of 
an anisotropic two-dimensional Ising model on a square lattice in a square geometry with periodic 
boundary conditions. Spins are coupled between nearest neighboring sites and between next-nearest 
neighboring sites along one of the lattice diagonals. We find that U* depends only on the asymptotic 
critical long-distance features of the anisotropy, irrespective of its realization through ferromagnetic 
or antiferromagnetic next-nearest neighbor couplings. We modify an earlier renormalization-group 
calculation to obtain a quantitative description of the anisotropy dependence of U* . Our results 
support our recent claim towards the validity of universality for critical phenomena in the presence 
of a weak anisotropy. 
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Universality is a key concept in the theory of criti- 
cal phenomena for bulk and confined systems, for re- 
views see, e.g., [1-3]. Within a given bulk universality 
class, critical exponents, certain critical amplitude ratios, 
and the critical behavior of thermodynamic functions are 
identical and independent of the specific microscopic re- 
alization. For instance, the bulk correlation length in the 
asymptotic critical domain, i.e., for asymptotically small 
positive or negative t = (T — T c )/T c , is, for isotropic 
systems, described by 

£ = £±,o|*r, T^T C , (1) 

where T c is the bulk critical temperature, v > is a 
universal critical exponent, £±fi are nonuniversal critical 
amplitudes, and = £+.o/£-,o is a universal critical 
amplitude ratio. 

Besides critical exponents, the thermodynamic func- 
tions describing the asymptotic critical domain are uni- 
versal if one allows for adjusting of only two amplitudes. 
Consider, for example, for small t and small field h con- 
jugate to the order parameter, the part /b. s of the bulk 
free energy density that becomes singular at the criti- 
cal point t, h = 0. /b, s is asymptotically described by 
universal scaling functions W ± for t ^ according to [4] 

/3/b,.(*, h) = A 1 \t\ 2 - a W ± (A 2 h\t\~ A ), (2) 

with f3 = l/(ksT), with universal critical exponents a 
and A, and where A\ and A 2 are nonuniversal ampli- 
tudes, i.e., they differ between different systems within 
the universality class under consideration. 

Universality also extends to situations, where the sys- 
tem under consideration is confined on a length scale L 
that is large compared to all microscopic length scales, 
such as lattice constants. If the singular part of the free 
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energy density of the system under consideration exhibits 
scaling, its asymptotic critical form may be written as [4] 

/3f s (t, h, L) = L~ d T{C x tl}l v , C 2 hL A /"), (3) 

where the nonuniversal constants C\ and C 2 are univer- 
sally related to the constants A\ and A 2 in (2). For given 
shape and boundary conditions, the function T is univer- 
sal. 

Equations (1)— (3) were stated for isotropic critical sys- 
tems, i.e., where the bulk asymptotic near-critical corre- 
lation length is isotropic. Here we are interested in the 
fate of universality in the presence of a weak anisotropy, 
i.e., where the amplitudes £±. n depend on the direction, 
but v does not (for strong anisotropics, where even v 
depends on the direction, see, e.g., [5] and references 
therein). Representatives of spatially anisotropic systems 
in the context of critical phenomena are, e.g., magnetic 
materials, alloys, superconductors [6], and solids with 
structural phase transitions [7, 8]. 

Most of the literature on critical phenomena has fo- 
cused on the isotropic case. In Rcfs. [9-18], investiga- 
tions were carried out on a number of specific anisotropic 
systems, but no general picture of universality for weakly 
anisotropic systems was suggested. The more recent pub- 
lications [19-23] have attempted to clarify the situation 
from a more general point of view. Refs. [19-21] con- 
cluded that the dependence of some bulk amplitude ra- 
tios, scaling functions, and the critical Binder cumulant 
on the parameters describing the anisotropy indicates a 
violation of universality. In Ref. [22] it was suggested that 
universality should be defined only after relating the sys- 
tem under consideration to an isotropic system by means 
of a shear transformation (see [9] for such a suggestion in 
the context of the two-dimensional Ising model). 

In contrast, the present author suggested in [23] that 
quantities that are universal in the isotropic case remain 
universal in the presence of a weak anisotropy, if their 
list of arguments is augmented by the parameters de- 
scribing the asymptotic critical long-distance features of 



2 



the anisotropy. Universality then implies that the quan- 
tity under consideration should exhibit no dependence on 
the particular microscopic realization of these features. 
Rather, it should depend universally on the anisotropy, 
i.e., in a way that is identical for all members of the re- 
spective bulk universality class. Explicit support for this 
claim was provided in [23] through results for free energy 
scaling functions of the two-dimensional Ising model on 
infinite strips. These examples, however, exhibited the 
simplification that a shear transformation exactly related 
the anisotropic systems to their isotropic counterparts 
with identical boundary conditions and geometry. Ac- 
cordingly, the scaling functions for the anisotropic case 
could be expressed in terms of the isotropic scaling func- 
tions, and the consequences of the anisotropy could be 
interpreted to be of mere geometric nature. In this work 
we consider a situation where this simplification is ab- 
sent. 

Consider weakly anisotropic critical bulk systems in d 
dimensions. Their asymptotic long-distance correlations 
are described by a correlation lengths ellipsoid that may 
be represented by a symmetric positive definite matrix 
S. This is diagonalized by a rotation matrix R, so that 

RER-^diagtf?, &...,$). (4) 

The £, are the asymptotic near-critical bulk correlation 
lengths along the principal axes of the ellipsoid. The 
scale-free matrix [23] 

A = (detH)- 1 ^ (5) 

is normalized to have det A = 1 and describes the shape 
and orientation of the ellipsoid. Since v is unique by as- 
sumption, A is independent of t within the asymptotic 
critical domain. It may be parameterized by d—1 correla- 
tion length ratios and d(d—l)/2 rotation angles. Accord- 
ing to [23] , universal quantities may receive an additional 
dependence on A, but are independent of its particular 
microscopic realization. This implies that the functions 
W ± (x) and T{y, z) on the right hand sides of Eqs. (2) 
and (3), respectively, have to be replaced by universal 
functions W ± (x, A) and T(y,z,A), which describe the 
isotropic situation as a special case A = 1. 

The use of a normalized matrix A to describe the 
anisotropy of a model was suggested in [19-21]. However, 
while A from (5) is defined through the physical corre- 
lation lengths in the asymptotic critical domain, explicit 
versions of A in [19-21] were obtained by expanding the 
Hamiltonian under consideration in small wave numbers 
k through order k 2 . While, for standard (p 4 field theory, 
these definitions coincide due to an exact mapping be- 
tween anisotropic and isotropic bulk Hamiltonians [19], 
this procedure in general does not, for lattice models, 
lead to the same A as defined through the physical cor- 
relation lengths. This will be demonstrated explicitly for 
the two-dimensional Ising model below. 

For d = 2-dimensional models, the ellipsoid reduces to 
an ellipse, whose shape may be characterized by the ratio 
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FIG. 1: (a) Bulk near-critical correlation lengths ellipse for 
two-dimensional models with an angle of inclination 9 and 
largest and smallest correlation lengths £> and £<, respec- 
tively. For the two-dimensional Ising model, the "1" and "2" 
directions are chosen as the m and n principal directions, re- 
spectively, of the square lattice, (b) Ising model square lattice 
with nearest-neighbor couplings J and next-nearest neighbor 
couplings J d on one of the diagonals. 

T = °f the smallest and largest correlation lengths 

(thus < r < 1) and an angle 9 describing its orientation 
(we select the convention — 7r/2 < 9 < n/2). If 6* is 
chosen to be the inclination of the direction of the largest 
correlation length with respect to the "1" direction, sec 
Fig. 1(a), the explicit form of A for d = 2 is [23] 

t _ / r sin 2 9 + cos 2 1 ,(r^ 1 -r)s'm29 \ . . 

§(r _1 -r) sin 29 r cos 2 + r" 1 sin 2 9) ' W 

Therefore, universal quantities may receive an additional 
dependence on the variables r and 9 as compared to the 
isotropic case. 

Consider a two-dimensional Ising model in a square 
L x L geometry on a square lattice with lattice constant 
a and ferromagnetic couplings J > of neighboring spins 
and couplings Jd of next-nearest neighboring spins in the 
direction of only one of the diagonals, see Fig. 1(b). With 
L = Na, the Hamiltonian of this model reads 

N 

H ^ ^ Sm,n [j(Sm+l,n "f" 5 m,n+l) ~t~ Jd s m+l,7i+l] > 

m,n— 1 

(7) 

with identifications of coordinates N + 1 and 1 in both 
principal lattice directions (i.e., periodic boundary con- 
ditions) and where s TOj „ = ±1 is the Ising spin at the site 
(m,n). We only consider couplings Jd > — J that allow 
for a ferromagnetic bulk phase transition at a tempera- 
ture T c > 0, given by [24, 25] 

sinh 2 (2/3 c J) + 2sinh(2/3 c J)sinh(2/? c Jd) = 1. (8) 

We choose 9 to be the inclination of the direction of the 
largest correlation length with respect to one of the lat- 
tice axes. Using results for the general triangular lattice 
[9, 23], we obtain, at the critical point, 

r = [sinh(2^ c J)] ±1 , = T7T /4 for J d ^ 0. (9) 
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Due to the symmetry of the L x L geometry, no depen- 
dence of universal critical quantities on the sign of 9 in (9) 
is possible. Thus, compared to the isotropic case, univer- 
sal quantities may receive an additional dependence only 
on r. Note that, for general r < 1, a shear transforma- 
tion to an isotropic system causes the periodicity of the 
boundary conditions to be no longer in mutually perpen- 
dicular directions and therefore no transformation to an 
isotropic system with identical geometry and boundary 
conditions to the original system is possible. Thus the 
consequences of the anisotropy cannot be interpreted to 
be of mere geometric nature. 

Here we test anisotropic universality for the critical 
Binder cumulant. The Binder cumulant [26] is a measure 
of the order parameter distribution. It may be used to 
locate the phase transition for a given model from the 
intersection of the cumulant for different system sizes or 
to determine the critical exponent v and therefore the 
universality class. For our system, the Binder cumulant 
is defined by 
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FIG. 2: (Color online) Transfer-matrix results from [30] and 
Monte Carlo results from [17, 18] for the Binder cumulant U* 
as a function of ln(l/r) [31]. Solid line: d=3 renormalization- 
group estimate for the relative change of U* for d—2 away 
from its isotropic value at r = 1. 



C/(T,i)^l-(M 4 )/(3(M 2 ) 2 ), 
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with the magnetization per spin M = N 2 J2 m n=\ a m,n 
and its second and fourth moments (M 2 ) and (M 4 ), re- 
spectively. We are interested in the critical thermody- 
namic limit U* = limi^oo U(T C , L). For given univer- 
sality class, geometric system shape (taken to the ther- 
modynamic limit), and boundary conditions, U* is be- 
lieved to be universal, i.e., to have a unique value in 
the isotropic case, independent of the underlying lattice 
[27, 28] and the spin value [29]. According to our con- 
siderations above, we expect U* to become a universal 
function U* (r) for the anisotropic case discussed here, 
i.e., if we switch on Jd. 

The most precise related results were obtained in [30] 
by means of a transfer-matrix method and imply 

U*(l) = 0.6106901(5), J d = 0, (11a) 
U*(l/V3) = 0.61182774(2), J d = J. (lib) 

Using Monte Carlo techniques, precision results have also 
been obtained for other values of Jd, notably for both 
,/ d > [17] and J d < [18]. All these results were dis- 
played in Fig. 4 of [18] as a function of s = 1/[1 + (J/Jd)]. 
Using the relation of r and s provided by the definition 
of s together with Eqs. (8) and (9), we have plotted all 
available data points from Refs. [17, 18, 30] as a function 
of ln(l/r) in Fig. 2. Note how all data points appear 
to be described by a single function U*(r), irrespective 
of the sign of Jd, i.e., for microscopic realizations of r 
through both ferromagnetic and antiferromagnetic cou- 
plings Jd. This is the central observation of this work, 
giving strong support to the conjectured universality of 
weakly anisotropic critical systems. 

We conclude by discussing the use of the renormal- 
ization group (RG) to describe the anisotropy depen- 
dence of U*. In [21], the RG in three dimensions was 



used to describe the anisotropy dependence of U* in 
two dimensions. A two-dimensional anisotropy matrix 
A.2(s) = (1 — s 2 ) -1 / 2 (I D with s from above was in- 
corporated into a block-diagonal 3x3 matrix A 3 (s) to 
estimate the s dependence of U* . The result provided a 
good approximation to the available data from [17, 30] 
for Jd > 0. However, later Monte-Carlo data for Jd < 
were not in agreement with the RG results, unless Jd is 
small [18]. 

The reason for this discrepancy is easy to identify. The 
matrix A.2{s) does not provide the correct description of 
the asymptotic shape of the correlation lengths ellipse in 
contrast to 



A 2 (r) = 



1 



r -+ 



r-f 



(12) 



obtained by setting 9 = = F 7r /4 in (6) for Jd ^ 0. While 
the deviation is small for positive Jd, it becomes large 
for negative Jd, unless Jd is small. Using A.2(r) instead, 
a natural choice for the appropriate three-dimensional 
matrix is 



A 3 (r) 




(13) 



which localizes the anisotropy within the first two di- 
mensions. Plugging our choice of A 3 (r) into the RG 
result for U^ G (r) = £7(0, A), provided in Eq. (7.8) of 
[21], and following Ref. [18] in normalizing the resulting 
function U^ G (r) by plotting U^ G (r)U*(l)/U£ G (l) with 
U*(l) from (11a), we obtain the solid line in Fig. 2, valid 
for both signs in (12) and thus for both ferromagnetic 
and antiferromagnetic J d . We conclude that the three- 
dimensional RG calculation of [21] results in a good de- 
scription of the relative anisotropy dependence of the 
two-dimensional critical Binder cumulant, provided the 
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anisotropy is described by the appropriate matrix A2(r) 
from (12) and A.2(r) is naturally embeddded into A^(r) 
as in (13). 

In summary, we considered the critical Binder cumu- 
lant U* of a two-dimensional Ising model on a square 
lattice in a square geometry with periodic boundary con- 
ditions and nearest-neighbor couplings. Additional cou- 
plings on one of the lattice diagonals allowed us to adjust 
a variable diagonal anisotropy of the correlation lengths 
in the asymptotic critical domain. We found that U* 
depends universally on its argument r, defined as the ra- 
tio of the smallest and largest bulk correlation length, 
independently of the particular realization of r through 
ferromagnetic or antifcrromagnctic diagonal couplings. 
Observing the true physical asymptotic long-distance be- 
havior of the anisotropy, represented by the asymptotic 
critical correlation lengths ellipse, we repeated Dohm's 



renormalization group calculation and found good quan- 
titative agreement of the relative anisotropy dependence 
of U* with precision transfer matrix and Monte Carlo 
data. Our results support the validity of universality for 
critical phenomena in the presence of a weak anisotropy. 
The ideas presented here and in [23] may be used to an- 
alyze other quantities in weakly anisotropic systems, for 
example those (besides the Binder cumulant), for which 
nonuniversality was claimed in [19-21]. 
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